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Analysis of full-genome sequences was previously used to trace the origin and transmission 
pathways of foot-and-mouth disease virus (FMDV) outbreaks in the UK in 2001 and 2007. 
Interpretation of these data was sometimes at variance with conventional epidemiological tracing, 
and was also used to predict the presence of undisclosed infected premises that were later 
discovered during serological surveillance. Here we report the genome changes associated with 
sequential passage of a highly BHK-21 -cell-adapted (heparan sulphate-binding) strain of FMDV 
arising from two independent transmission chains in cattle. In vivo virus replication rapidly 
selected for a wild-type variant with an amino acid substitution at VP3^®. Full-genome sequence 
analysis clearly demonstrated sequence divergence during parallel passage. The genetic diversity 
generated over the course of infection and the rate at which these changes became fixed and 
were transmitted between cattle occurred at a rate sufficient to enable reliable tracing of 
transmission pathways at the level of the individual animal. However, tracing of transmission 
pathways was only clear when sequences from epithelial lesions were compared. Sequences 
derived from oesophageal-pharyngeal scrapings were problematic to interpret, with a varying 
number of ambiguities suggestive of a more diverse virus population. These findings will help to 
correctly interpret full-genome sequence analyses to resolve transmission pathways within future 
FMDV epidemics. 



INTRODUCTION 

Foot-and-mouth disease (FMD) is a highly contagious, 
acute viral disease of cloven-hoofed animals, characterized 
by fever, loss of appetite, depression, lameness and the 
appearance of vesicles on the feet and in, or around, the 
mouth. Recovery generally occurs within 8-15 days; 
however, ruminants can carry FMD virus (FMDV) in the 
oropharynx for years following the resolution of acute 
infection (Alexandersen et al, 2002). The virus can spread 
extremely rapidly (partly as a consequence of the small 
amount of virus that can initiate infection, the large 
amount of virus excreted by affected animals and the 
multiple routes of infection) by direct or indirect contact 
with infected animals or their products or by long-distance 
airborne transmission (Alexandersen et al., 2003). FMDV 
has the potential to cause enormous economic losses and is 
the single most important constraint to international trade 
in livestock and animal products (Perry & Rich, 2007). 

The GenBank/EMBL/DDBJ accession numbers for the FMDV 
sequences included in this study are JX570638-JX570655. 

A supplementary table is available with the online version of this paper. 



FMDV belongs to the genus Aphthovirus (family 
Picornaviridae) . The FMDV virion consists of an icosahed- 
ral shell composed of 60 copies each of the four structural 
proteins, VP1-VP4 (Rueckert & Wimmer, 1984). The 
capsid proteins surround a positive-sense ssRNA genome 
approximately 8.3 kb in size that replicates via a negative- 
strand intermediate (Belsham, 2005). FMDV, in common 
with other RNA viruses, displays very high mutation rates 
during replication, corresponding to the poor fidelity of 
the RNA polymerases (Domingo et al, 2001). The 
mutation rate, together with other factors, for example 
the genomic architecture, replication speed and recomb- 
ination, determines the rate at which a virus population 
evolves (Duffy et al., 2008). Commonly cited mutation 
rates for RNA viruses lie in the range 10~^-10~^ mutations 
nt~' per genomic replication (Drake & Holland, 1999). 
FMDV RNA polymerase is no exception, with current 
estimates suggesting that, on average, one nucelotide 
sequence change occurs during each round of virus 
replication (Haydon et al., 2001). In the field, the high 
mutation rate, fast replication cycle and large size of the 
affected population results in the rapid evolution of FMDV 
(Domingo et al, 2003). Molecular epidemiology studies to 
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monitor the evolution of the relatively short nucleotide 
sequences encoding the capsid VPl protein provide 
valuable insight into the emergence of new strains and 
serotypes worldwide (Ayelet et al, 2009; Valarcher et al, 
2009). However, full-genome sequencing of FMDV pro- 
vides greater resolution when reconstructing transmission 
pathways, as demonstrated in the UK epidemic in 2001 and 
in 2007 (Cottam et al, 2006, 2008b). These reports are 
based on consensus genomes sequenced directly from 
epithelial lesions, blood or oesophageal-pharyngeal scrap- 
ings (probang samples) collected fi-om individual animals 
on infected premises (IPs). The 2001 post-outbreak 
investigation retrospectively resolved inter-farm transmis- 
sion of FMDV (Cottam et al, 2006, 2008a; Konig et al, 
2009). In 2007, full-genome sequencing was used in 'real 
time' to identify the initial outbreak virus and to connect 
the first and second temporal clusters of IPs (Cottam et al, 
2008b). These data were invaluable for informing control 
policies, they predicted the existence of undetected 
intermediate premises that were subsequently identified, 
and demonstrated that full-genome sequence analysis has 
the required resolution for inter-farm transmission tracing 
(Cottam et al, 2008b; Schley et al, 2008). The initial 
outbreak virus was identified as Oi/BFS 1860/UK/67, a 
strain propagated extensively in cell culture for vaccine 
production (Cottam et al, 2008b). Field isolates of FMDV 
use integrin heterodimeric (a//?) glycoproteins as cellular 
receptors through an interaction mediated by the Arg-Gly- 
Asp (RGD) motif within VPl (Jackson et al, 2000). Serial 
passage of FMDV in cell culture may select for viruses that 
also bind to heparan sulphate through acquisition of 
positively charged amino acid residues at VP2'^'* and/or 
VPS^" (Fry et al, 1999; Jackson et al, 1996; Sa-Carvalho 
et al, 1997). An additional change from a negatively 
charged amino acid at VPS*'" to a neutral residue often 
occurs; however, this change is not considered essential for 
heparan sulphate binding (Sa-Carvalho et al, 1997). These 
heparan sulphate-binding phenotypes are considered to be 
attenuated in host species (Sa-Carvalho et al, 1997). 
Herein we report the genome changes associated with 
reversion to in vivo replication for a highly cell-adapted 
strain of FMDV during sequential, parallel passage in a 
controlled study in cattle. Despite having undergone 
extensive passage in cell culture, in vivo virus replication 
in the natural host rapidly selected for a wild-type variant 
with an amino acid substitution at VP3^*. We demonstrate 
that FMDV transmission can be reliably traced at the level 
of the individual animal, with clear divergence of sequences 
during sequential, parallel passage. 

RESULTS 



FMDV transmission and disease presentation 

Pyrexia, nasal discharge and large vesicles on the feet and 
within the oral cavity, secondary to the inoculation sites, 
were present within 24 h of intradermolingual challenge 



(Figs 1 and 2). All calf-to-calf challenges resulted in clinical 
FMD, with epithelial vesicles detected in the oral cavity, 
tongue and on all four feet of all animals. Vesicular fluid or 
epithelium from ruptured vesicles was collected immedi- 
ately upon detection, or within 24 h if detected before 
calf-to-calf challenge in order to prevent artificial envir- 
onmental contamination. The samples selected for sequen- 
cing (as described in Methods) are shown in Fig. 2. 

Maximum-parsimony analysis and transmission 
patliways 

Three phylogenetic trees were constructed, including: (i) all 
the sequences determined (not shown); (ii) only the 
sequences from viruses isolated from animals in the A 
group (Fig. 2b); and (iii) the sequences from viruses isolated 
from animals in the B group (Fig. 2d). The transmission 
events based on sequences isolated from epithelial lesions 
were clear and unambiguous, with one to four nucleotide 
changes for each calf-to-calf transmission: IAH2 (a tissue- 
culture-adapted FMDV with a heparan sulphate-binding 
phenotype) to Al (1 nt), Al to A2 (2 nt), A3 to A5 (1 nt), 
IAH2 to Bl (4 nt), Bl to B2 (2 nt), B2 to B3 (4 nt) and B2 to 
B4 (1 nt). However, some events were complicated by 
apparent back-mutations: A2.6D.V to A3.5D.V (Fig. 2b) and 
B3.3D.V to B5.9D.V (Fig. 2d). Nevertheless, this may 
indicate that these particular virus isolates were not part of 
the direct line of transmission and the apparent back- 
mutations actually existed in another site in the same 
animal. The combined tree clearly indicated that the virus 
infecting calf A4 was derived from calf B4 (three nucelotide 
differences as opposed to 10 if it had originated from donor 
calf A2). Calves A4 and B4 (indirect contacts with calves A2 
and B2, respectively) were housed together following 
challenge (Fig. 1) and FMD lesions were recorded for calf 
B4 3 days before calf A4, providing a window during the 
period of direct contact for virus transmission. 

Sequences derived from oesophageal-pharyngeal scrapings 
were problematic to interpret, with a varying number of 
ambiguities in all sequenced samples and fi-equently 
diverging from the main branch of transmitted virus (Fig. 
2). Interestingly, the viral sequences derived from calf A2 and 
B2 oesophageal-pharyngeal scrapings at 32 days post- 
challenge (p.c.) were phylo genetically more closely related 
to the inoculum than to viral sequences derived from 
epithelial lesions from both animals at 6 days p.c. (seven to 
nine changes compared with 10 to 13 changes, respectively). 
This suggests that virus may evolve independently in different 
sites in the same animal and this may be more obvious in the 
oesophageal-pharyngeal area, where the virus can establish 
infection and persist for long periods. 

Distribution and analysis of evolutionary changes 
at the consensus level 

A total of 49 nucleotide substitutions at 47 sites along the 
genome were detected in 18 consensus sequences (Table 1). 
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Fig. 1. FMDV infection and calf-to-calf transmission. Phase 1 : two calves (A1 and B1) were challenged by subepidermolingual 
injection and housed in the same box for 24 h. Phase 2: each inoculated calf was moved into a separate box to challenge a 
naTve calf (A2 or B2) by direct-contact exposure. Phase 3: after 5 days of direct-contact exposure, vesicular lesions were 
detected in the oral cavity, on the tongue, snout or mouth area of the contact-exposed calves. Inoculated calves A1 and B1 
were removed from the study. Calves A2 and 82 were moved to the portioned boxes, each containing two naive calves, one for 
direct-contact challenge (A3 or 83) and a second for indirect-contact challenge (A4 or 84). Phase 4: following a 24 h 
challenge, the two indirect-contact challenged calves (A4 and 84) were removed from the portioned boxes and housed 
together for 14 days. The direct-contact challenged calves (A3 and 83) were moved to separate boxes containing a naive calf 
(A5 or 85) for an additional direct-contact challenge for 14 days. Calves A2 and 82 were moved into a single box and 
maintained up to 32 days p.c. (phase 5). 



The genetic variation observed between viruses was mainly 
due to synonymous point substitutions, the majority of 
which were transitions. The non-synonymous substitutions 
resulted in 10 amino acid substitutions (Table 2). The 
transition : transversion ratio for the A transmission chain, 
excluding sequence A4.12D.V, was 6.67. The ratio for the B 
transmission chain, which included A4.12D.V, was 3.67, 
demonstrating a strong bias towards transition substitutions. 
The 3' UTR and 3B1 and 3B2 coding regions contained the 
greatest number of substitutions per site (0.056, 0.029 and 
0.028, respectively; Table 1). There were no changes detected 
in the VP4, 2A or 3B3 coding regions. One to three 
ambiguities at the consensus level were detected in all of the 
sequence traces derived from oesophageal-pharyngeal scrap- 
ing samples. In comparison, no ambiguities were detected in 
sequence traces derived from epithelial lesions. 

In vivo passage of FMDV IAH2 rapidly selected for wild- 
type variants, with substitution at residue VP3^^ from Arg 
to Cys detected at the consensus level in both of the 
inoculated calves at 2 days p.c. (Tables 1 and 2). This initial 
selection was conserved at the consensus level through all 
four parallel calf-to-calf transmissions. However, an 
oesophageal-pharyngeal scraping from donor animal A2 
at 32 days p.c. (A2.32D.P, Fig. 2) contained a variant virus, 
with a substitution at VP3^* from an Arg to a His residue 
(Tables 1 and 2). No substitutions were detected at 
positively charged Lys at VP2^^* or neutral Gly at VP3^". 
Besides the change at nucleotide position 2771 (resulting in 
a Cys at VP3^^), only a substitution of C to U at nucleotide 
position 170, within the 5' UTR S-fragment, was present in 
both group A and B (Table 1). In addition, there was one 
disparate substitution at position 8157 within the 3' UTR 



in both groups. All other changes occurred at different 
nucleotide positions, demonstrating clear sequence diver- 
gence between the A and B transmission chains. 

DISCUSSION 

Cattle passage of IAH2, a tissue-culture-adapted FMDV 
with a heparan sulphate-binding phenotype, caused clinical 
FMD within 24 h and rapidly selected for a wild-type 
variant. Reversion at the consensus level was detected 2 days 
p.c. in both inoculated calves by substitution at VP3^* from 
Arg to Cys. This initial selection was conserved at the 
consensus level through all four calf-to-calf transmissions on 
two parallel transmission chains. However, an oesophageal- 
pharyngeal scraping collected from one calf at 32 days p.c. 
contained a variant with a substitution at VP3^* from an Arg 
to a His residue. No substitutions were detected at VP3''°, a 
residue less likely to undergo reversion (Fry et al., 1999; Sa- 
Carvalho et al, 1997), or at VP2'^^''. 

Previously reported next-generation sequence analysis of 
IAH2 revealed only low-frequency Cys and His nucleotide 
coding substitutions at the VP3^* coding region (Table SI, 
available in JGV Online) (Wright et al., 2011). Single 
nucleotide polymorphisms (SNPs) at frequencies between 
0.17 and 0.06% at this site are indistinguishable from 
background, consistent with errors inherent in RNA 
replication prior to onset of any selective pressure (Wright 
et al, 2011). No high-frequency SNPs were detected in either 
the A or the B group transmission chains following parallel 
in vivo passage. The observed substitutions were either 
random or the result of mutations which rapidly selected 
against heparan sulphate binding. Clearly, extensive culture 
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Fig. 2. FMD progression and phylogenetic trees, (a, c) Tabular representations of the experiment (Fig. 1), highlighting samples 
selected for sequencing. Each column of the table under phases 1-5 represents a study day. FMDV infection and disease 
presentation are colour-coded for each animal: red shading represents initial day of challenge, light-blue shading represents 
incubation period before clinical FMD, orange shading represents period of clinical FMD and green shading represents day of 
termination. Samples are coded AX or BX.YD.V or P: AX or BX, animal number for group A or B, respectively; YD, days post- 
challenge (p.c); V, vesicular epithelium or vesicular fluid; P, probang sample (oesophageal-pharyngeal scraping sample), (b, d) 
Phylogenetic trees constructed from the viral sequences derived from the A or B groups, respectively. Samples are colour- 
coded as described above, with lines between circles representing a single nucleotide change. Samples that contained 
consensus-sequence ambiguities are shown as larger circles (labelled with the actual number of ambiguities that were present), 
which are linked to other samples via the closest sequence contained within the population. Sample A4.1 2D.V is represented by 
a grey circle and dashed line (from A1 .2D.V) in (b), as the transmitted virus originated from animal B4 and is part of the B 
transmission group (d). 



passage and subsequent adaptation to bind heparan sulphate 
did not completely attenuate IAH2 in vivo. It was previously 
demonstrated that a chimeric FMDV Oi Campos/Ai2 cDNA 
clone, with Arg at VPS^*^, was highly attenuated for cattle at 
lower doses (Sa-CarvaUio et al, 1997). However, reinfection 
with a higher dose resulted in clinical FMD, and recovered 
virus had reverted back to wild type with an uncharged Cys 



at VP3 or a negatively charged Glu at VP2 (Sa-Carvalho 
et al, 1997). A His at VP3^'' was also detected following in 
vivo passage (Sa-Carvalho et al, 1997). The involvement of 
Arg at VP3^* in viral attenuation for cattle has also been 
confirmed using Oi/Campos/Bra/58 cDNA clones (Borca 
et al, 2012). Similar studies in swine also resulted in 
delayed onset of FMD, with recovered virus exhibiting 
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Table 1. Distribution of evolutionary clianges (as shown by nucleotide substitutions in viral proteins) at the consensus level 



The sequences derived from oesophageal-pharyngeal samples are in bold. These sequences contained ambiguities and frequently diverged from the main branch of the transmitted virus (Fig. 2). 
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*Samples are coded AX or BX.YD.V or P: AX or BX, animal number for group A or B, respectively; YD, days p.c; V, vesicular epithelium or vesicular fluid; P, probang sample (oesophageal- 
pharyngeal scraping sample). 

tThis row indicates the number of nucleotide substitutions per site. 
tThis row indicates the genome position for each nucleotide substitution. 

§The virus infecting calf A4 (A4.12D.V) was determined to have originated from calf B4 (B4.9D.V); therefore, A4.12D.V is part of the B transmission group. 
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Table 2. Distribution of amino acid clianges in viral proteins following sequential parallel passage in cattle 



The sequences derived from oesophai 


;eal-pharyngeal samples are in bold. These sequences contained ambiguities and frequently diverged from the 


main branch of the transmitted 


virus (Fig. 2). 






Viral protein 


affected by amino acid substitution 










T pflHpr VP? 




VP3 


2C 


3A 


3C 


PP338* 


PP356 PP364 PP659 


PP894 


PP961 PP1570 


PP1887 


PP2107 


P4t 


P22 P30 P39 


P56 


P123 


P129 


P128 


P124 


1104$ 


1158 1181 2066 


2771 


2772 2972 


4799 


5750 


6410 


Sequence§ 
















IAH2 


T 


L K A 


R 


D T 
K 1 


A 

A 


T 


V 


A1.2D.V 






C 










A2.2D.P 


I 




C 










A2.4D.P 






c 










A2.6D.V 






c 










A2.6D.P 






c 










A2.32D.P 




E/KII 




H 








A3.5D.V 






c 










A4.12D.V? 






c 




s 






A5.9D.V 






c 










dL.ZU. V 






c 










B2.2D.P 




P/LII 


C/RII 










B2.4D.P 




P/LII T 


c 




s 






B2.6D.V 






c 




s 






B2.6D.P 




P 


c 


V 


S/AII 






B2.32D.P 




P/LII 


c 


V 




A 




B3.3D.V 






c 




s 




I 


B4.9D.V 






c 




s 






B5.9D.V 






c 




s 




I 


Charge# 


None 


None + to — None 


+ to N 


+ to N None 


None 


None 


None 



""This row indicates the polyprotein position of the amino acid substitution. 
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§Samples are coded as explained in Table 1. 

I Ambiguous nucleotide substitutions. 

?The virus infecting calf A4 (A4.12D.V) was determined to have originated from calf B4 (B4.9D.V); therefore, A4.12D.V is part of the B transmission group. 
#Change in charge as a result of non-synonymous nucleotide substitution: None, no change; -I- , positive; — , negative; N, neutral. 
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reversion to Cys or His at VP3 (Borca et al, 2012). 
Interestingly, these authors also demonstrated that in vitro 
heating of the viruses at 41 C caused a pronounced loss of 
infectivity of VP3^'' Arg phenotype viruses compared with 
VP3^* His viruses (Borca et al, 2012), which may play a role 
in selective pressure in vivo. Here we show how rapidly a 
single nucleotide mutation can be selected during in vivo 
passage, and that reversion at VP3*° and VP2^^* are not 
required for virulence and therefore presumably also for wild- 
type receptor (integrin) usage. 

The UK 2007 FMD outbreak was initiated by escape of OJ 
BFS 1860/UK/67 from the Pirbright complex, leading to 
two IPs in August and to a second cluster of IPs 6 weeks 
later and approximately 20 km apart (Schley et al, 2008). 
Three closely related viruses, lAHl, IAH2 and MAH, used at 
the Pirbright site were identified as possible source viruses 
(Cottam et al, 2008b). The heparan sulphate-binding strains 
MAH and IAH2 are closely related to virus isolated from the 
first IP (IPlb, six and seven substitutions, respectively). In 
comparison, lAHl virus, which lacks residues associated 
with heparan sulphate binding, differs by 12 substitutions 
(Cottam et al, 2008b). MAH can be distinguished from 
IAH2 by one substitution at position 1181 (A to G, a non- 
synonymous change at amino acid residue 2 of the Leader-b 
polypeptide) which was also present in all 1 1 reported field 
virus sequences (Cottam et al, 2008b). We have shown that 
IAH2 virus can cause clinical FMD in cattle, with rapid 
reversion to wild-type (non-heparan sulphate binding) Cys 
or His at VP3^*. Sequential passage resulted in accumulation 
of nucleotide substitutions but did not provide strong 
selective pressure for reversion at VP3*'', consistent with all 
but one of the 2007 field viruses (Cottam et al, 2008b). 
These data are consistent with evidence that a cell-culture- 
adapted virus was responsible for the 2007 outbreak. 
Significantly, we demonstrate clear viral sequence divergence 
during sequential passage through parallel groups of cattle, 
even in the face of a selective pressure for receptor usage and 
in vivo replication. Only changes at nt 2771 (VP3^*) and 
change C to U at nt 170 were present in both the A and the B 
transmission chains. The evident divergence of sequences 
and accumulation of nucleotide substitutions provides 
support for the interpretation of molecular epidemiological 
analysis during 2007 that linked the two outbreak clusters, 
contrary to available field epidemiological evidence (Schley 
et al, 2008). Cottam et al (2008b) identified 36 nucleotide 
differences between IAH2 and 11 UK-2007 field isolates. 
Interestingly, we observed seven of these substitutions 
following in vivo passage of IAH2, although these changes 
were not present as SNPs above 1 % in the inoculum 
(Wright et al, 2011). Only a single change at position 170 
was common to both the A and the B transmission chains. 
Taken together, these data may suggest that position 170 was 
under strong selective pressure during cattle passage. 
However, the significance of this change is not clear, as it 
is predicted to convert a G-C base pair of a conserved stem 
structure into a G-U base pair within the 5' UTR S-fragment 
(data not shown). 



We have shown that the genetic diversity generated over 
the course of infection, and the rate at which these changes 
became fixed and were transmitted between calves, 
occurred at a rate sufficient to enable reliable tracing of 
transmission pathways at the level of the individual animal. 
However, transmission events were only clear and 
unambiguous when based on sequences derived from 
epithelial lesions. Sequences from oesophageal-pharyngeal 
scrapings contained ambiguities, complicating their inter- 
pretation, and often diverged from the main branch of 
transmitted virus. The viral diversity in oesophageal- 
pharyngeal scrapings may be the result of continuous, 
low-level virus replication at this site, compared with the 
transient period of replication, prior to sampling, at sites of 
vesicular lesions. Alternatively, the diversity observed may 
simply reflect the lower levels of virus present in these 
samples, resulting in amplification errors. Further studies, 
utilizing next-generation genome sequencing technology, 
are required to understand the intra-host population 
dynamics of FMDV and to identify the source of viruses 
for onward transmission to other animals. 

Surprisingly, the one to four nucleotide changes that we 
detected during controlled inter-animal transmission were 
comparable to the number of changes detected between 
inter-farm transmissions during the UK FMD epidemic in 
2001 and in 2007 (Cottam et al, 2008b). During the 2001 
analysis, 23 consensus sequences were recovered directly 
from epithelial lesions acquired from 21 IPs over approxi- 
mately 7 months, and the nucleotide changes were estimated 
to accrue at a mean rate of 1.5 substitutions per farm 
infection (Cottam etal, 2006) . Additional analysis by Cottam 
and colleagues of 22 consensus sequences from 15 IPs over 
approximately 3 months indicated a mean of 4.3 substitu- 
tions upon farm-to-farm transfer (Cottam et al, 2008a). 
During the 'real-time' analysis in 2007, 11 consensus 
sequences were recovered directly from epithelial lesions, 
oesophageal-pharyngeal scrapings or blood samples acquired 
from 1 1 IPs over approximately 2 months, and the consensus 
sequences between farms differed by one to five substitutions 
(Cottam et al, 2008b). It is conceivable that the small number 
of nucleotide changes that we detected during controlled 
inter-calf transmission is in part a consequence of the 
relatively short period that FMDV-infected cattle are 
infectious (mean 1.7 days) (Charleston et al, 2011). The 
short infectious period, combined with the small amount of 
virus required to initiate infection and multiple routes by 
which the virus can spread, limits the time available for the 
virus to accumulate consensus-level substitutions before it is 
transmitted and establishes infection in the subsequent host. 
In addition, the rapid replication rate combined with a rapid 
host immune response reduces the susceptibility to super- 
infection with closely related strains (Juleff etal, 2009). These 
inter-animal transmission characteristics of FMDV could 
also account for the small number of changes detected 
between IPs in the UK during 2001 and 2007. The small 
number of substitutions between IPs suggests that the contact 
network between farms was very high, and there were limited 
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cycles of intra-farm animal-to-animal transmission prior to 
disease detection and sampling. More in-depth analysis of 
multiple samples fi^om each IP is required to understand the 
intra- and inter-farm population dynamics of FMD during 
the 2007 outbreak. In particular, the possible impact of 
bottlenecks on transmission dynamics is not clear, for 
example intra-farm transmission cycles initiated by indirect 
contact and subsequent low quantities of source virus. 

Our results suggest that, following transmission of virus to 
a subsequent herd, continual intra-farm transmission cycles 
on the IP or ongoing virus replication in a single host during 
persistent infection, as detected in oesophageal-pharyngeal 
scrapings, can generate discrete viral populations that are 
genetically divergent from the transmitted virus. Analyses of 
sequences derived from these samples are difficult to 
interpret and can be misleading, with implications for 
correct interpretation and resolution of transmission path- 
ways. Therefore an understanding of intra-farm viral 
dynamics is particularly important when there is a broad 
spectrum of lesion ages on an IP or when oesophageal- 
pharyngeal scrapings are the only sample type available for 
tracing. These data demonstrate that fuU-genome sequen- 
cing can be used to resolve FMDV transmission events at the 
level of the individual animal and highlight the importance 
of the sample type and the description of the clinical 
situation for accurate interpretation of transmission 
events. 



METHODS 

Inoculum. The FMDV isolate Oi/BFS 1860/UK/67 (IAH2) (GenBank 
accession no. EU448369) used for challenge was originally derived 
from bovine tongue epithelium received at the Food and Agriculture 
Organization of the United Nations (FAO) World Reference 
Laboratory for FMD at Pirbright in 1967 from a farm near 
Wrexham, England. After extensive passage in baby hamster kidney 
(BHK-21) cells, it was found to have adapted to utilize heparan 
sulphate as a cellular receptor (Cottam et al, 2008b). 

Experimental design. The experiment was designed for sequential, 
parallel passage of FMDV through two groups, A and B, each 
containing five calves (Fig. 1). The experiment consisted of five phases 
carried out in high-containment large animal isolation boxes at the 
Pirbright Laboratory. Two calves (Al and Bl) were challenged with 
FMDV by subepidermolingual injection of 10^ '' TCID50 (phase 1) as 
described previously (Henderson, 1952). Each inoculated calf was 
moved 24 h p.c. into a separate box to challenge a naive calf (A2 or 
B2) by direct-contact exposure (phase 2). Vesicular lesions were 
detected in the oral cavity of calves A2 and B2 5 days p.c. Inoculated 
calves Al and Bl were then removed from the study and each 
contact-exposed calf (A2 and B2) was moved to a separate portioned 
box to act as a donor to transmit FMDV over 24 h by direct or 
indirect-contact exposure (phase 3). Each portioned box contained 
two naive calves, one for direct-contact challenge (A3 or B3) and a 
second separated by a 30 cm wide, double-walled wooden-and-mesh 
partition for indirect-contact challenge (A4 or B4). During the first 
6 h, the ventilation in the portioned boxes was turned off, raising the 
relative humidity to >99% (Charleston et al, 2011). In addition, a 
wall-mounted fan was used to direct air from the donors to the 
indirect-contact calves over the 24 h period. Following challenge, the 
two indirect-contact challenged calves (A4 and B4) were removed 



from the portioned boxes and housed together for 14 days (phase 4). 
The direct-contact challenged calves (A3 and B3) were moved to 
separate boxes containing a naive calf (A5 or B5) for an additional 
direct-contact challenge for 14 days (phase 4). Calves A2 and B2 were 
moved into a single box and maintained up to 32 days p.c. (phase 5). 
This was necessary since British Home Office rules precluded the 
housing of lone animals for extended periods and space was limited. 
Calves were examined daily during phases 1^ for signs of clinical 
FMD. Animal experimentation was approved by the Institute for 
Animal Health (lAH) Ethical Review Board under the authority of a 
Home Office project licence in accordance with the Home Office 
Guidance on the Operation of the Animals (Scientific Procedures) 
Act of 1986 and associated guidelines. 

Sample collection and coding for sequencing. Samples collected 
and sequenced were coded AX or BX.YD.V or P: AX or BX was the 
animal number for group A and B, respectively; YD was the number of 
days p.c; V was for epithelium or vesicular fluid collected from FMD 
epithelial lesions and P was for oesophageal-pharyngeal scrapings 
(probang samples). Vesicular fluid was collected using sterile syringes, 
diluted to a 10% solution with M25-phosphate buffer (35 mM 
Na2HP04 . 2H2O; 5.7 mM KH2PO4; pH 7.6; made in-house) and 
stored at — 70 °C until used for RNA extraction. Vesicular epithelium 
was collected using sterile forceps and tweezers, immersed in cryotubes 
(Sarstedt) containing 50 % (v/v) glycerol in M25-phosphate buffer and 
stored at —70 °C until being processed. Approximately 1.5 g vesicular 
epithelium was ground by using a pestle and mortar and resuspended 
to a 10% solution with M25-phosphate buffer. The solution was then 
centrifiiged for 10 min at 3500 g at room temperature, and the 
supernatant was removed for RNA extraction. Oesophageal-pharyn- 
geal scrapings collected using probang cups (http://www.wrlfmd.org/ 
find_diagnosis/probang.pdf) were stored at —70 °C until used for RNA 
extraction. 

Vesicular epithelium was collected 2 days p.c. from the hoof of the 
left pelvic limb of calf Al (A1.2D.V) and from the hoof of the right 
pelvic limb of calf Bl (B1.2D.V). Oesophageal-pharyngeal scrapings 
were collected from calves A2 and B2 on days 2, 4, 6 and 32 p.c. 
(A2.2D.P, A2.4D.P, A2.6D.P, A2.32D.P and B2.2D.P, B2.4D.P, 
B2.6D.P, B2.32D.P) and vesicular epithelium was collected 6 days 
p.c. from the hoof of the left pelvic limb of calves A2 and B2 (A2.6D.V 
and B2.6D.V). Vesicular epithelium was collected 5 days p.c. from the 
hoof of the right thoracic limb of calf A3 (A3.5D.V) and vesicular fluid 
was collected 3 days p.c. from the hoof of the right thoracic Umb of calf 
B3 (B3.3D.V). Vesicular fluid was collected 12 days p.c. or 9 days p.c. 
from the hoof of the left thoracic hmb of caff A4 (A4.12D.V) and B4 
(B4.9D.V), respectively. Vesicular fluid was collected 9 days p.c. from 
the hoof of the right thoracic limb of calf A5 (A5.9D.V) and from the 
hoof of the left thoracic limb of caff B5 (B5.9D.V). 

Genome amplification, sequencing and sequence analysis. 

Viral RNA was extracted directly from clinical samples using either an 
RNeasy Mini kit (Qiagen) or TRIzol reagent (Invitrogen) as described 
previously (Reid et al, 1998). Reverse transcription and DNA 
amplification were performed as described previously (Cottam et al, 
2008b), but with a pre-amplification cDNA purification step using an 
fflustra GFX PGR DNA and Gel Band Purification kit (GE Healthcare 
UK). Negative controls were used for each step, including control 
reactions for each of the 24 primer sets used for full-genome 
amplification. The quantity of viral RNA used was not standardized 
between the samples that were sequenced. 

Amplified PGR products were separated by gel electrophoresis, 
stained with ethidium bromide and visualized under UV light. After 
direct purification, cycle sequencing was carried out as described 
previously using either a Beckman DTCS kit (Beckman Coulter) on a 
Beckman Coulter CEQ 800 sequencer (Cottam et al, 2008b) or a 
BigDye Terminator v3.1 Cycle Sequencing kit (Applied Biosystems) 
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on an ABI PRISM 3730 analyser. The SeqMan Pro program from 
the Lasergene v8.0.2 package (dnastar) was used for assembly 
and proofreading of sequence trace files. Alignment and manipulation 
of sequences were performed using BioEdit v7.0.5.3 (Hall, 1999) 
and DnaSP v4.10.3 (Rozas et al, 2003) software. Maximum- 
parsimony analysis of the nucleotide sequences was performed using 
TCS vl.21 (Clement et al, 2000) as described previously (Cottam et al., 
2008b). Where ambiguities were present, multiple sequences that 
contained the different substitutions, but no ambiguities, were 
constructed; using these sequences, TCS generated complex multi- 
linked trees. The shortest routes through these trees (which did not 
involve back-mutations) were chosen to construct the final trans- 
mission trees. 
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